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A variational formulation is given for a theory of gravity coupled to a massive vector in four dimensions, 
with Asymptotically Lifshitz boundary conditions on the fields. For theories with critical exponent z — 2 
we obtain a well-defined variational principle by explicitly constructing two actions with local boundary 
counterterms. As part of our analysis we obtain solutions of these theories on a neighborhood of spatial 
infinity, study the asymptotic symmetries, and consider different definitions of the boundary stress tensor 
^^ . and associated charges. A constraint on the boundary data for the fields figures prominently in one of our 

JH , formulations, and in that case the only suitable definition of the boundary stress tensor is due to Hollands, 

Ishibashi, and Marolf. Their definition naturally emerges from our requirement of finiteness of the action 
under Hamilton- Jacobi variations of the fields. A second, more general variational principle also allows the 
Brown- York definition of a boundary stress tensor. 
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1 Introduction 

An interesting extension of AdS/CFT is the study of models that arise in condensed matter physics, particularly 
those describing quantum critical systems [1, 2, 3, 4]. Such models, which are difficult to study using traditional 
methods based on weakly interacting quasiparticles and broken symmetry, may exhibit non-relativistic scaling 



symmetries of the form 

t^XU x^Xx (1.1) 

with dynamical critical exponent z ^ 1. These symmetries, present at quantum critical points, provide a strong 
kinematic connection to some versions of the AdS/CFT correspondence. The transformations (1.1), known as 
Lifshitz scahng, are encoded in the asymptotic symmetry group of the conjectured gravitational dual theory. 

An early manifestation of this idea is based on four-dimensional gravity with a vector and a 2-form that 
interact via a topological coupling [5]. However an equivalent (and simpler, for our purposes) formulation is 
obtained by integrating out the 2-form [6]. The result is gravity coupled to a massive vector, described by the 
Lagrangian 



^=— V=5|"^-2A--F^^F^,- — ^'^AJ . (1.2) 



2^2 V .y 4 ^'^ 2 

This Lagrangian yields the following equations of motion for the metric and vector 

4d 1„ 4d_ a^ i^/^Z7A77. 1^ 77AK77. \ : '^ f A A ^„ ^A 



'Rf^u - 2 5/.. ""R = -A 5^. + - [Ff,^ F,x - - gt^u F^^Fx^j + — ^Af,A, - - g^, ^M^ j (1.3) 

^V^'F^, = n?A,. (1.4) 

With the appropriate boundary conditions, the equations of motion admit solutions that realize the non- 
relativistic scaling symmetries. The simplest solution of this kind is 

ds'= i^S dr' - i"X dt' + {-X dx^ (1.5) 



A,dx^= i^y f-^dt. (1.6) 

The critical exponent z and the length scale i that characterize the solution are related to the parameters of 
the theory by 

/2z . z'^ + z + A 



For z = 1, the vector field vanishes and the metric (1.5) describes the Poincare patch of AdS4. Other solu- 
tions with the same boundary conditions describe localized excitations of a massive vector (i.e., excitations 
with compact support) on an asymptotically AdS4 spacetime. The asymptotic symmetry group of the properly 
formulated theory is the three-dimensional conformal group SO(3,2). However, when z ^ 1 the symmetries 
change. The asymptotic symmetries of these 'Lifshitz' solutions [7, 8] include the non-relativistic scaling trans- 
formations (1.1) (noted previously in other contexts [9]), translations, and spatial rotations, but not boosts. 
The case z = 2 is of particular interest, because certain models of strongly correlated electrons are invariant 
under the transformations t — ;► A^ t, x — >■ A x. 

As in AdS/CFT, the on-shell action is an important tool in studying the properties of the dual Condensed 
Matter models. However, simply integrating the Lagrangian (1.2) over spacetime and adding the Gibbons- 
Hawking- York term [10, 11] at spatial infinity 
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I=^Jd''x./^rR-2A--F'"'F„,,,-^A''A„,\+^l d'xV-hK (II 




does not give an appropriate action. An immediate problem with this action is that it exhibits divergences 
and other undesirable properties when evaluated on solutions asymptotic to (1.5)-(1.6). This is because - as 
is often the case when one makes the jump from (1.2) to (1.8) - the EOM are not actually stationary points 
of the proposed action. The variation of this action contains surface terms that do not vanish unless the field 
variations satisfy more restrictive fall-off conditions than the solutions themselves. 

This is a common problem when formulating gravitational theories on non-compact spacetimes, where a 
clearly defined boundary value problem is not automatically equivalent to a well-defined variational principle. 
For example, a variational principle for asymptotically flat gravity requires a new boundary term in addition to 
the usual Gibbons-Hawking- York term. Otherwise the action is not stationary under generic 1/r deformations 
of the Schwarzschild solution. This was originally understood from the Hamiltonian point of view [12], and 
more recently investigated for the Lagrangian description of the theory [13, 14, 15, 16, 17]. 

In this paper we give a proper variational formulation for a theory of gravity coupled to a massive vec- 
tor with "Asymptotically Lifshitz" solutions. We start by determining a set of boundary conditions for the 
fields that generalizes the Lifshitz asymptotics and includes spacetimes with curved spatial sections. We then 
determine the proper definition of the variational principle for the theory by supplementing the action (1.8) 
with an appropriate set of surface terms (boundary counterterms) , a procedure that has come to be known as 
Holographic Renormalization [18, 19, 20, 21, 22, 13, 23]. This extends the results of [24], which investigated 
the variational formulation of theories with solutions that asymptote to (1.5)-(1.6). 

There are two significant assumptions in our approach. First, a completely general analysis of the equations 
of motion is quite difficult, so we make a few simplifying assumptions - guided by the form of known solutions - 
about the asymptotic behavior of the fields. While less general, this renders the analysis tractable, and highlights 
key features that must be present in a full treatment of the problem. Second, we focus on gravitational theories 
that provide a dual description of CM models with dynamical exponent z = 2. Thus, we do not claim to have a 
completely general definition of an Asymptotically Lifshitz spacetime (even for the specific case z = 2). While 
our results are quite broad, a more complete analysis of asymptotically Lifshitz boundary conditions (including 
time-dependence of the fields) is left for a future work. 

The outline of our paper is as follows. In section 2 we formulate the equations of motion in a manner 
that lends itself to the analysis of spacetimes with anisotropic dependence on r, give boundary conditions 
that generalize the asymptotics of (1.5)-(1.6) to solutions with curved spatial sections, and then construct 
solutions that satisfy these boundary conditions locally on a neighborhood of spatial infinity. We also describe 
the asymptotic symmetries of the theory. In section 3 we demonstrate that the action (1.8) is not suitable 
for the variational formulation of the theory, and then construct an action with the requisite properties. The 
construction is based on two requirements: that 5/ = on-shell for field variations having the same asymptotics 
as solutions to the equations of motion, and that bl is finite on-shell for variations of boundary data that respect 
the kinematic constraints of the theory. Finiteness of the on-shell action follows from these more fundamental 
conditions. We then extend this construction and obtain an action such that 5/ = on-shell for any field 
variations that preserve the boundary conditions of the theory. In section 4 we use the procedure introduced by 
Hollands, Ishibashi, and Marolf (HIM) in [25] to construct the "improved" boundary stress tensor and conserved 
charges for the theory. We also show how one of the actions introduced in section 3 allows the more familiar 
Brown- York (BY) stress tensor to be defined for the theory. The charges constructed from the BY stress tensor 



are also conserved, but only as a result of some fairly restrictive assumptions in our boundary conditions. The 
HIM charges are related in a suggestive way to the BY charges and a component of the boundary current dual 
to the massive vector field. We explicitly demonstrate how to make use of our results in section 5, using Lifshitz 
topological black holes [7] as an example, and close our paper with a discussion in section 6. 

Concerning notation, our calculations involve two ADM-type decompositions: a 4 — > 3 + 1 decomposition 
for studying the evolution of the fields along a spacelike coordinate r, followed by a 3 —^ 2 + 1 decomposition 
that splits the fields into parts with distinct asymptotic behavior. This requires notation for quantities on four-, 
three-, and two-dimensional spaces, which we summarize here. The spacetime is a four-dimensional manifold 
Ai with metric g^i, and coordinates x^. A superscript '4' is prepended to spacetime curvatures and covariant 
derivatives, as in (1.3)-(1.4). Three-dimensional surfaces M C Ai have coordinates x"" and metric hab- The 
intrinsic curvatures on these surfaces are denoted Rat and R, the covariant derivative is Va, and the extrinsic 
curvature is Kab- Finally, two-dimensional spatial surfaces T, G M have coordinates x* and metric aij. The 
Ricci tensor on S is TZij, the covariant derivative is Di, and the extrinsic curvature is 9ij. 

2 Asymptotically Lifshitz Spacetimes 

In this section we generalize the asymptotics of the solution (1.5)-(1.6) to spacetimes with curved spatial sections. 
We begin by expressing the equations of motion in a form suitable for studying these solutions. This is a two- 
step process, consisting of a 3 + 1 split that emphasizes the evolution of the fields along the spacelike coordinate 
r, followed by a 2 + 1 split that separates the fields into components with distinct asymptotic dependence on r. 
The resulting equations are used to identify a simple set of 'asymptotically Lifshitz' boundary conditions, and 
then solved in a neighborhood of spatial infinity to determine the asymptotic behavior of the fields. 

2.1 Decompositions of Fields and Equations of Motion 

For the 3 + 1 decomposition, we work with coordinates x^ = (r, x") adapted to a foliation {Mr} of the spacetime 
by constant r surfaces. It is convenient to partially fix the coordinate gauge so that the metric in the asymptotic 
region r ^ i takes the form 
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g^u dx^'dx" = ^dr^ + habix", r) dx^dx^ , (2.1) 



with hab{x'^,r) the metric on a surface Mr- Tensors are projected normal or parallel to Mr by appropriate 
contractions with the unit vector n'^ = {r/l)5^r or the projector Pa^ = dx'^/dx"'. For the vector A^ this gives 
a normal component An = "ny Ay and a three-vector Aa = Pa^{Afj_), while the components of its field strength 
F^iu are 

Pa^Pb'iF^u) = daAb - dbAa := Fab (2.2) 

Pa'in^F^,) = £nAa - daA^ := Ba , (2.3) 

with £n the Lie derivative along n'^. The projections of the equations of motion are carried out in the same 
manner, which gives five equations involving the fields hab, Aa, An, and their derivatives. The first three 



equations, obtained from the Einstein equations (1.3), are 



1 (i^2 _ i^afe^^, j _^R= _A+jB''Ba - ^ F"''F,b + ^ V " ^ ^"^'^ ^^'^^ 



where 



VfoK"^ - V"K = — ^n ^"^ + - F^^Bfo (2.5) 

Gab + hab£nK - £nKab+ 2 Ka'Kbc " i^Ti^Tab + ^ /iafe ( i^' + K'^^K^d ) (2.6) 

= -/l,b A + ^ Fa'^Fkc - ^ Kb F'^'^F.d + \ BaB, - ^ hah BcB'' 

9 9 9 

-f^afe = -^ £nhab (2.7) 

is the extrinsic curvature of M, K = /i''^ i^af, is its trace, and Gab is the three-dimensional Einstein tensor. The 
two remaining equations are the projections of the Proca equations (1.4) 

-VaB'' = m^An (2.8) 

VFab + £nBb - 2 BaK\ + KB^ = m^Ak . (2.9) 

In analogy with the usual ADM decomposition, (2.4), (2.5), and (2.8) are "constraint equations", while (2.6), 
(2.9), and the definitions (2.3, 2.7) of Ba and Kab are "evolution equations". These 3 + 1 equations do not rely 
on any assumptions other than the choice of coordinate gauge. 

The next step is a 2 + 1 split of the three-dimensional coordinates into a time coordinate t and spatial 
coordinates x^ , which gives a fohation {Sj} of each M^ by spatial surfaces of constant t. For a general choice 
of coordinates x" = {t,x^) the 3-metric takes the form 

habdx'^dx^ = -a^ dt^ + aij {dx' + /3' dt){dx^ + /3^ dt) , (2.10) 

where a is the lapse function, /3* is the shift vector, and o",-,- is the spatial 2-metric on a surface Ej. The 3- vector 
Aa admits a similar split into a temporal component </> and a spatial 2-vector Ai 

Aadx" = (pdt + Ai [dx' + P'dt^ . (2.11) 

In principle, the 3 + 1 equations (2.4) - (2.9) should now be projected normal and parallel to S^ to yield a set of 
'2 + 1 + 1' equations of motion. However, implementing this for a completely general set of fields and solving the 
resulting equations is quite complicated, so we restrict our attention to fields whose asymptotic form is given 

by 

habdx^dx'' = -a^{x^,r)dt^ +a,j{x^,r)dx'dx^ (2.12) 

A^dx^ = (/){x'',r)dt . (2.13) 



That is, for r ^ I the fields are assumed to be independent of t, and certain components - the shift vector /?*, 
the scalar An, and the 2-vector Ai ~ may be ignored. These restrictions let us avoid many of the complications 
of the general 2 + 1 split, while still accommodating the original solution (1.5)-(1.6), the topological black hole 
spacetimes of [7], and other solutions. They need not hold on the full spacetime Ai, and our results should 
apply equally well to more general field configurations when the deviations from (2.12)-(2.13) are suppressed in 
the asymptotic region. 

The details of the 2 + 1 decomposition of the equations of motion are discussed in appendix B. The evolution 
equations associated with the definitions of Kab and Ba split cleanly into temporal and spatial components 

Ktt = -a £nOi Kti = Kij = - £ncrij (2-14) 

Bt = £n^ B, = . (2.15) 

The various projections of (2.4)-(2.9) yield a set of nine equations, four of which are automatically satisfied by 
fields of the form (2.12)-(2.13). The remaining equations are 

1 {K\f - \ K'^Kij + - £naK\ -\n + -D^a = -K--^ (£„0)2 + -1^ Ac^-I^V + ^ '/'' (2.16) 

2 la z a 4a^ 4a^ 4a^ 



^ "'^^ ^ ^ ^ ^^ ' -£n4> (k\ - - £na) - — </> = . (2.18) 

a. \ n ' a. 



Di(-D'(t>]+-£n'(l) + 

a ^ a a ^ a ^ a 
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DjK'^ - D'K^j = dH- £na) - K'^ - DjU + — ^ ^ndD'a ^ £n(t>D'(j) (2.19) 



Ih" 



a — DiDja j + aij — £n a + CFij£nK'^rn — £nKij + 2 Ki^Kjm — K'^m Kij 



+ -£„«( GijK'^m - K,j) + \ a,, {K'^mf + \ Oij K^'^Krnk (2.20) 

11 1 m^ 

where indices are lowered and raised using the 2-metric aij and its inverse. While the 3 + 1 equations were 
completely general, these 2 + 1 + 1 equations of motion are only valid when the fields have the restricted form 
(2.12)-(2.13). 

2.2 Boundary Conditions and Asymptotic Behavior of Solutions 

Boundary conditions that generalize the asymptotics of the solution (1.5)-(1.6) can be determined from a 
straightforward analysis of the 2 + 1 + 1 equations of motion. The leading behavior of the fields for r ^ £ is 
assumed to take the form 

ai,{x\r) ~ (0'^S'(^') «(x^r) ~ {^-^\''\x') <^(x^r) ~ {^jj r\x') , (2.21) 



where crf^ , a*''\ and </>*"' are smooth functions on S^ that comprise the boundary data for the fields. Consistency 
with the equat: 
the constraint 



with the equations of motion does not place any restrictions on 0"-°', but the functions q*"' and (/)'"' must satisfy 



(/>t°' = ± aW W ^ ^^ ^\ (2.22) 

Invariance of the theory under t ^ —t means that the choice of sign in this equation is not important; we will 
always assume the positive sign. Thus, boundary data for the theory consists of a spatial 2-metric a^^^ and a 
single scalar function that we define as a'"' = e^^. The asymptotic behavior (2.21) with smooth boundary data 
that satisfies (2.22) will be referred to as asymptotically Lifshitz (AL) boundary conditions. 

Solutions of the equations of motion with AL boundary conditions can be constructed locally on a neigh- 
borhood of spatial infinity for any choice of a^j and x- There is no guarantee that these solutions also exist 
globally, but we will only require properties of solutions that are present in the asymptotics. First we perform 
a Taylor expansion of each field, starting from the leading behavior (2.21) 



2-n 



E 7 4"'(-') « = E 7 «'"'(-') ^ = E 7 ^'"'(-^). (2.23) 



(Tij 

The equations of motion are then solved order-by-order in powers of i/r, which determines the coefficients in 
the expansions (as local functions of the boundary data) up to some finite order that depends on z. Beyond 
this point the coefficients also depend on dynamical aspects of the fields that are not fixed by the boundary 
conditions. 

For the rest of this paper we will focus exclusively on theories with critical exponent z = 2. In this case the 
expansions are qualitatively similar to an asymptotically AdSs spacetime [21, 22], with vanishing coefficients 
for the n = 1 and n = 3 terms in (2.23), and the coefficients of the n = 4 terms only partially determined by 
the equations of motion. The asymptotic expansions of the fields out to this point are^ 

{-^ + {^-^ + {^-^ + --) (2-24) 

a= (y%^-(l+(^)'a^^) + (^)'a<^) + ...) (2.25) 

where we have pulled out a convenient overall factor of e^-^ in a and (/>, and use '. . .' to represent higher-order 
terms (n > 4) that will not be needed for our analysis. The n = 2 and n = 4 terms will be referred to as 
next-to-leading order (NLO) and next-to-next-to-leading order (NNLO), respectively. The equations of motion 



^For AL boundary conditions with z — 2 the expansions do not contain log{r/£) terms. We have not ruled out such terms for 
different boundary conditions, or critical exponents z > 2. 



determine the NLO terms as two-derivative combinations of the boundary data 

<'' = -y <°' ^^"' + ^' 2).2),x + l^ 'Q^xT^a - ^ <^S' 2)'X - y 4'^ 2)'=X23,X (2.27) 

a(=^> = :^Da23*X (2.28) 

</>(^> = y 7e<«' + ^1)2;^ + ^ ©,X25^;^ , (2.29) 

where !Dj is the two-dimensional covariant derivative compatible with ai~- and TV-^^ is its scalar curvature. Thus, 
the NLO terms are fixed by the kinematicas of the theory. But at NNLO the fields are only partially determined 
by the boundary data. A single scalar degree of freedom appears, which we take to be the NNLO term 0*'*' in 
the expansion of the vector field. The NNLO term in the expansion of the lapse, a*'*^ , can be written in terms 
of c/)*''^ and four-derivative combinations of the boundary data as 



5 " 320 ^ ^ 80 80 ^80 



«(4) = _ i ^(4) _ L^ ^TiW f -L 'D^'jz'o) _ ^ ^.j^(o) ^t^ + l-7^(«) V'xV.x (2.30) 



p4 pA Of A f4 of^- 



-j^g A A i ^Q A A JA ^g A "^ 40 "^ 40 

Both the trace and divergence of the NNLO term in the spatial metric are fixed by the equations of motion. 
Like a^'*' , the trace a^^^ = a^^ a^f depends on the scalar degree of freedom 
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5 ^ 320 ^ ^ 160 80 ^10 ^ ^ ^ ' 

/4 0/4 f4 ni o/)4 

- 1^ D2;, T)\ - ^ 1)2;^ l^^xO^a - ^ 2)^X 2)'2)a + 1^ 2)22)2;, + :^ ( V\V,X? 

+ ^ 2)*x2)^X 2).2),X + ^- W^x 'D^'DJX ■ 
The divergence is determined by a complicated equation of the form 

V^alf + 2 4^^ V^x + V, = , (2.32) 

where Vj is a combination of LO terms, NLO terms, and their derivatives. Solutions of this equation include a 
homogenous part annihilated by the first two terms, which may be written in terms of a transverse tensor Yij 



^(4)hom ^ e-2xy.^.. The full result for a^;^ is 
a§ = e-'^Y,, - i a^ cP^^^ + ^ ag' (7^'") )^ - ^ ( 2)aI^,7^^"^ + D.7^<»)D,x) - ^ 2),2),7^(°) (2.33) 

+ ^ ag> D27^(°) - ^ a^ 7^(°) D^x - y 'Da'D.x ^^'x - y ( 2?a 2?'2),X + 2),X D^^a ) 

/)4 »4 7/)4 /)4 

- i^ 2)2a3,a^,.X + ^ 4?' 2)\ 23.7^'°) - {|^ ag> 7^('') V^x ^kX - y D.2),X D^x2?.X 



— a,^. D X"^ VkX + ^'y^, D 2) X+^Q ^y ,^ a^.a; . « ^y 



+ ^ -.^?^ ^'X 1>'V,x + ^ cxg' D^D^x + ^ ^g' ( 2)^X 2).x )' + V ^^ 2)^D'x l^.^a 



The trace of the inhomogenous part of this solution agrees with (2.31), so in addition to being transverse the 
tensor Yij must also be traceless. We show in appendix C that symmetric, transverse-traceless tensors on S^ 
built from four-derivative combinations of the boundary data vanish identically, which confirms that Yij is 
independent of the boundary data. 

The expressions for the fields out to NNLO are quite complicated, so it is useful to look at a special case 
where the results simplify. In section 5 we consider the topological black hole solutions [7], which are examples 
of AL fields where the function x ~ the boundary data for a and (/> - is constant. For simplicity we will take 
X = 0, in which case the coefficients in the asymptotic expansions (2.24)-(2.26) are 

(2) _ _^ (0) Tp(o) (4) _ Y ._- o-'"' ci)'*' + — a^°^ f7^'°' )^ -—VV •7^("^ + — ct'"' D^T?'") (2 341 

% ~ 8 ^ij '^ ^ij ~ -^»J ^^ij ^ + 640 *^ ^ ' 32 ^^ 160 ^-^ ^ ^ ^ 

a(^)=0 aW = -ic/,w__L(7^(o))2_L 1)27^(0) (2.35) 

5 320 ^ ^ 80 ^ ^ 

^(2) ^L 7^(0). (2.36) 

8 

As with the more general solution, the degrees of freedom appearing at NNLO are a scalar 0'"' and a transverse- 
traceless tensor Yij. 

To summarize, AL boundary conditions for the restricted set of fields (2.12)-(2.13) are given by the leading 
asymptotic behavior (2.21), with boundary data consisting of a scalar function and a metric on some two- 
dimensional spatial surface S^. Solutions of the equations of motion can be constructed locally on a neighborhood 
of spatial infinity, and take the form (2.24) - (2.26) for theories with critical exponent z = 2. The boundary 
data at LO fixes the NLO terms according to (2.27) - (2.29), but the NNLO terms are only determined up 
to a scalar degree of freedom and the transverse-traceless part of the spatial metric. At higher orders in the 
expansion, the coefficients depend on this dynamical information in addition to the boundary data. 

2.3 Asymptotic Symmetries 

The asymptotic symmetries of a theory are the diffeomorphisms that act on the fields at spatial infinity, while 
preserving the boundary conditions. For a theory with boundary conditions that admit the Lifshitz solution 
(1.5)-(1.6), these include the non-relativistic scaling transformations (1.1). For other choices of AL boundary 
conditions the Lifshitz scaling transformations are typically broken, though other symmetries may exist. 

Under an infinitesimal coordinate transformation x^ ^ x^ — S,^{x) the metric and vector transform as 

^i9tj.u = £^gij.u S^Afj_ = £^A^ . (2.37) 

In order to preserve the coordinate gauge (2.1), these diffeomorphisms must take the general form 

C = rjix'^) e = a)(x'=) - / dr^ h'^' 9,7 , (2.38) 

where 7 and ^?pj may depend on the three-dimensional coordinates x'^, but not on r. Requiring that the 
diffeomorphisms also respect the asymptotic form of the fields (2.12)-(2.13) further restricts the coordinate- 
dependence of the components of ^^, and we find 

C = r^{x') e = CUit) e = eUix') + Q' ^Xx'^) + (^y el,,{x') + ... . (2.39) 



The NLO and NNLO terms in ^* are obtained using the asymptotic expansion of a*-^ in the integral in (2.38), 
which gives 



e 



j.) = ^I^'7 6J,, = -^42),7. (2.40) 



Thus, diffeomorphisms that preserve both the coordinate gauge and the asymptotic form of the fields consist 
of a rescaling of r that depends on the spatial coordinates x^ , a t-dependent reparameterization of time, and a 
diffeomorphism on Sf that includes contributions that are sub-leading in r. 

The action of the diffeomorphisms on the fields can be worked out from the transformations (2.37). We are 
primarily concerned with the response of the fields at LO and NLO in the asymptotic expansion, since these 
terms are completely fixed by boundary conditions and kinematics. Under the diffeomorphism (2.39), the LO 
and NLO terms in the expansion (2.24) for the spatial metric transform as 

<%^l°' = 27^i°' + i^e(o,ag' (2.41) 

5^a^ = £,^a^ + £,^af:^ . (2.42) 

Now consider the fields a and cj). Before the constraint (2.22) is applied, the asymptotic expansions take the 
form 

The terms in the expansion of the lapse transform according to 

5^a(°) = (27 + a* e[o)) «'"' + -fejo)"'"' (2.44) 

5^a^^^ =-27 a(^' + £,(Q, a*'' + £,^^.^ log a'^' , (2.45) 

and likewise for the terms in the expansion of (j) 

'.(°'= (2^ + dti%)<P''' + £e,,,ct>'~'' (2.46) 

(2) = _ 27a(^) + £,^^^0(^) + £,^^^ logcA*"' . (2.47) 

The overall factors of a*°' and 1;^*"' in (2.43) result in slightly non-standard transformations for the NLO terms. 
Notice that the transformations of the LO terms preserve the constraint on the boundary data (2.22), so we 



conditions 



can restrict our attention to a'°' and af^ when determining the diffeomorphisms that preserve the boundary 



Asymptotic symmetries act at spatial infinity but preserve the boundary conditions, which means that 
the transformations of the LO terms in the fields must vanish. In the case of the lapse, requiring ^^a*"' = 
immediately restricts the possible form of the time reparameterization. Since 7 and a'°' are functions of x , it 
follows that 9(^*0) cannot depend on t if (2.44) is to vanish, and the only possibilities are constant rescalings or 
constant translations 

^*^^ = -2At-(5t . (2.48) 

10 



Taking this into account, the conditions for the diffeomorphism to preserve the boundary data can be written 
as 

£,,„,a(°' = (2A + Dfce(;Ja(o) (2.49) 



£,^^a!^=a^V,e'^,,, (2.50) 



a'°' and al°\ asymptotic symmetries are associated with A and ^*qj that satisfy these equations. For the original 
Lifshitz solution (1.5)-(1.6) we have x = ^^'^ '^ij = ^iji ^^^ in that case the asymptotic symmetries include 



where we have used the trace of (2.41) to express 7 in terms of the divergence D^^^^gj. For given boundary data 

^j ■ , asynipLULic syninieirieB are associaieu wim a anu q,^^ +1-"+ „„^-,„<-,, j-v,„ 

olution (1.5)-(1.6) we he 
the Lifshitz scaling generated by 

^^ = Ar e* = -2At e = -^x'. (2.51) 

Other choices of boundary data may or may not allow non-trivial solutions of (2.49)-(2.50). Notice, however, 
that the asymptotic symmetries always include the constant time translation t ^ t + 5t. 

One application of our results is to AL solutions with a'°' = (j)'-°^ = 1 and spatial sections of constant (non- 
zero) scalar curvature. In these cases, the asymptotic symmetries do not include Lifshitz scaling transformations. 
The easiest way to demonstrate this is to consider the transformation of the NLO terms in the fields. Since 
the NLO terms are completely fixed by the kinematics, they should not change if the diffeomorphism preserves 
the boundary conditions. Using the x = solutions (2.34) - (2.36), the transformations of the fields at NLO 
become 

6^a^^^ = (2.52) 

'^«'/'''^ = ^(4A-2A)7^'°) (2.53) 

hc7^= -^4"'(6f„,5,-2A)7^'°), (2.54) 

where we have simplified some terms and canceled others using (^^a'"' = and S^a^^ = 0. If the scalar curvature 
T?.'"' is a non-zero constant, then the NLO terms in the vector and spatial metric only vanish if A = 0. The 
equations (2.49)-(2.50) might have non-trivial solutions in these cases, but they do not include the Lifshitz 
scaling transformations (1.1). 

3 A Variational Principle 

A variational principle identifies solutions of a theory as stationary points of an action. Generally speaking, 
there is some space of allowed field configurations, and the action must be stationary for any variation of the 
fields within the space. A basic requirement for a "well-defined" variational principle is that this space should 
include generic field configurations with the same asymptotic behavior as any physically reasonable solution of 
the theory [12]. If a proposed action is not stationary for arbitrary variations within this space, then it is not 
suitable for the variational formulation of the theory. Of course, one may have additional applications of the 
action in mind, which require that it have certain properties on an even larger space of field configurations. 
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The results of the last section can be used to make the basic requirement on the space of field configurations 
more precise for theories with AL boundary conditions and critical exponent z = 2. In that case, all solutions 
are identical at LO and NLO in the asymptotic expansion, but they may differ at NNLO. This means that 
an action - when evaluated on a solution of the equations of motion - must be stationary for independent 
variations of the fields at NNLO. As we pointed out in the introduction, the action (1.8) does not have this 
property. The response of that action to a small change in the fields is 

SI = ^ f d^'x./^ (£^'''6g^, + £^'6A^+^ f d^xV^iir'^'dhab + TT^SAa) , (3.1) 

where S'^'^ = and £^ = give the equations of motion (1.3)-(1.4), and the coefficients of the field variations 
in the surface integral are 

tt"'' = - (h"^ K - K"'''^ 7r'^ = -ln^F^V (3.2) 

The bulk integral in (3.1) vanishes for solutions of the equations of motion, but the surface integral does not. 
To see why this is the case, it is convenient to work with the fields a, aij, and </>. The surface integral is now 
given by 



61 



\ I d?x V^ (-2 a ^" 6a + tt'^ 6aii + vr* 6(1)) . (3.3) 

K^ JdM ^ ' 



^-° '^^ 'dM 



To make sense of this integral, which contains factors that either vanish or diverge at spatial infinity, the 
integrand should be evaluated on a regulating surface M^, with r ^ i. The asymptotic expansions for the 
fields can then be used to determine the r — > cxd limit. The field variations fall off as r~^ in this limit ^, since 
they behave like NNLO terms in the asymptotic expansions (2.24)-(2.26). But their coefficients - including the 
contribution from the volume factor - grow as r^. As a result, the surface integral makes a finite but non-zero 
contribution to 61 as Mj- is taken to spatial infinity^. We conclude that the action (1.8) is not stationary under 
the full class of variations required for a well-defined variational principle. 

3.1 A Minimal Action 

An action with the appropriate variational properties is obtained by adding new surface terms to (1.8). Since 
the action is a functional of fields that diverge at spatial infinity, we work on a compact region Mr '^ M 
bounded by a surface Mr of constant r ^ i. The action is then defined as the r — t- oo limit of a functional Ir 
on this cut-off spacetime. We assume that the new surface terms are at most quadratic in the fields and their 
derivatives, so Ir is given by 



Ir = ^ J^d^x ./^ rR - 2A - -F'^'^ F^, - — Af'A^] + -^ jd-'xV-hK (3.4) 




+ ^ / d^X V^ (co + Cl A'^Aa + C2R + C3 F'^^Fab 

K'^ J Mr ^ 



^This is not the case if we work with htt instead of a. The variation Shu ~ —2a 5a approaches a constant at spatial infinity. 
^Notice that the field variations themselves go to zero in this limit. A Dirichlet boundary value problem is not equivalent, in 
this case, to a well-defined variational principle. 
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Any calculation involving the action or its variation is performed using this functional, and completed by the 
r — )■ oo limit that takes M^ to spatial infinity. Other procedures for cutting off the spacetime may give different 
results, so our choice of limiting procedure should be thought of as part of the definition of the theory. 

The coefficients of the surface terms in (3.4) are fixed by demanding a specific response from the action 
for two different types of field variations. First, 51 must vanish on-shell for small, independent variations of 
the fields at NNLO. This guarantees that solutions of the equations of motion are stationary points of the 
action within the basic space of field configurations needed for a well-defined variational principle. Second, 
the action should have a finite response to field variations caused by small changes in the boundary data that 
respect the constraint (2.22). This insures that both the on-shell action and the conserved charges (obtained 
from a suitably defined boundary stress tensor) are finite. It also has the effect of enlarging the space of 
field configurations allowed by the variational principle. An action that meets these two requirements will be 
referred to as 'minimal'. It is worth pointing out that the choice of surface terms in (3.4) is not unique, and 
other combinations of surface terms can be used to obtain a minimal action. This issue will be discussed in 
more detail at the end of this section. 

We will now use the conditions on the variation of the action to determine the coefficients {q} in (3.4). 
Since the new surface terms are intrinsic to the regulating surface M^, they do not introduce additional bulk 
terms in 51. The on-shell variation of the action has the same basic form as before 



^^1-0 -^X/-^ 



(3.5) 



[tt'^'+P^') 5K, + {^''+p'')5Aa 
with the contributions from the new surface terms given by 

p"'' = i Co h"^ + ci Q h"^ A^A^ - A''A^\ + C2 Q /i"^ R - R""^) + cg Q h"^ F"^F^d - 2 F\ F^A (3.6) 

p"" = 2ciA''-Ac3VbF^'' . (3.7) 

For the variational principle to be well-defined, (3.5) should vanish for independent variations of the fields at 
NNLO. These can be written as 

5a^ 



'ij 



(^) H'^ ^"=(^) e'"''^"*" '^'^=(^) ^'"''^'^^'^ (3.8) 



with the factors of e^-^ in 5a and 5(j) included for convenience. Using these expressions in 51, we have 

2 / n\2 



51 



\ f d^x ^fH \ {n'^ + p'^) (-) 5a'if - 2 a (vr" + p") (-) e^^ 5a^'^ (3.9) 

K J Mr L V ^ / V '^ / 



^-° K" JM,. 



For r ^ £, the r~'^ behavior of each field variation is canceled by the r"^ growth of its coefficient, leaving three 
independent terms in the surface integral that are finite and non-zero as r — )■ oo. For the action to be stationary, 
the coefficients {q} must be tuned to cancel these finite contributions. Working out the asymptotic expansions 
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for each term in the integrand, the leading behavior is 

2 



V-h (vr'J + p'^ 
-2v^^a(7r"+p" 



-H-'+p'){^'' 






\2l^ 2 



(0) p2x 



'cjw e 



(0) p2x 



2 

Co , Ci 



Co Ci \ jj (4) 



'(0)""jj 



'criu; g 



+ - + -Ua- + ... 



2ci <5(/)(^' + ... , 



(3.10) 
(3.11) 
(3.12) 



with '. . .' indicating terms that vanish as r — > oo. Only the first two surface terms that were added to (3.4) 
make finite contributions to 51 for the field variations (3.8); the other terms fall off too rapidly at large r. This 
poses a potential problem, since there are two coefficients and three independent terms that must be addressed. 
However, all three terms are canceled by setting 



1 



Co 



Cl 



(3.13) 



Thus, requiring the on-shell action to be stationary for independent variations of the fields at NNLO determines 
the coefficients of two of the surface terms in the action. 

The coefficients of the remaining surface terms are fixed by demanding a finite response of the action to 
Hamilton- Jacobi variations of the fields; i.e., field variations associated with small changes of the boundary data 
that satisfy the constraints of the theory. These variations take the form 



5(Ti 



-A < + ... 



5a 



, r 
2 - 



=2X 



5x + ... 



2 [-\ e^>^5x + ... . 



(3.14) 



Notice that the variations of a and (p have the same leading behavior, since the kinematic constraint (2.22) 
M"' when z = 2. The '. . .' in each expression is a reminder that the NLO terms (2.28) - (2.27) 



requires a'°' 

will also change when the boundary data is varied. However, the contributions to 51 from the NLO terms in 
the field variations will vanish in the r — t- oo limit, so we may ignore them in this calculation. The change in 
the action due to a Hamilton-Jacobi variation of the fields is 



51 



1 



d^x 



-h 



TT*-' + p 



Wl 



50^^ + - 4 a (tt" + p«) + 2 (tt* + p 



e^^5x 



. (3.15) 



Working out the asymptotic expansion of the terms in the integrand reveals contributions proportional to 
r^, which diverge as r — >■ oo, and contributions independent of r, which are finite in that limit. The terms 
proportional to r"^ can be removed by tuning the coefficients of the remaining surface terms. First we consider 



the term in (3.15) proportional to Sa- : 



(0). 



-h {tt'^ + p'^) 



^-^ 



1 /r 

4 



3^^Vo<^ 



- 16c2 - 32 C3 - 21) {V'x-D'x - 4) '^''x'DkX) (3.16) 
+ ( - 8C2 - 2^) {V^V^x - 4) '^''^kX) 



+ {i-16c3)a%V^x1)kX 



5a. 



(0) 



14 



with '. . .' indicating the part that is finite as r — >■ oo. The r"^ terms are canceled by setting 



C2 



'4 



C3 



16 



(3.17) 



These values of C2 and C3 also cancel divergences in the r^ terms proportional to 6x, which are given by 

2 



-h 



4 a iir'' + p") + 2 (vr 



..'))G 



e^^6x 



(3.18) 



=2X 



/cr(0) 



(4 C2 + ^'°^ + (16 C3 - (D^25,,x + ^''x'DkX) 



Sx + 



Therfore, the on-shell response of the action to Hamilton-Jacobi variations of the fields is finite as r — )■ cxd if the 
coefficients of the surface terms are given by 

5 I £ i , , 

co = -^ ci = - C2 = -- C3 = -. (3.19) 

Although it is not obvious from the expression (3.18), the fact that ^a*"' and ^c/)*"' satisfy the constraint (2.22) 
on the boundary data is essential to canceling some of the r^ terms in the variation of the action. If we repeat 
this calculation with independent variations Ja'"' and ^c/)*"' , the r^ terms in 51 cannot be canceled for any choice 
of the coefficients {c,} in (3.4). This is examined in more detail in the next section. 

With the coefficients of the surface terms fixed, the action is 



2^2 



d^x- 



Mr 



-g{''R-2A 






'^A>'A, 



+ 



K? J Mr 



d^xV^K 



(3.20) 



+ ^ / d'^x 

K^ J Mr 



-h 



Ad A 



16 



-R + — F^'Fab 



In addition to having the required behavior under the different variations of the fields, this action is also finite 
as r —7- 00 for solutions of the equations of motion. The exact value of the action will depend on the particular 
solution, but the common asymptotics of the fields is enough to show that the on-shell action does not contain 
surface terms proportional to positive powers of r. Using the equations of motion and integration-by-parts, the 
on-shell action can be written 



1 
2^2 Jj^^ 



d'^x yf^ \2K + 



m 



A^'A,, 



+ -J d\ 

K'^ J Mr 



-h 



1 



n. F^'^ A, + K 



(3.21) 



+ 



K^ Jm, 



d^xV^ 



Aa A 



4^+16 



pabj, 



The surface terms are evaluated in the asymptotic region r ^ i, making potential divergences easy to identify. 
The bulk term, on the other hand, depends on the behavior of the fields on the whole spacetime. Luckily, 
divergences associated with the r — >■ 00 limit may be identified by replacing the integrand with its asymptotic 
behavior, performing the integral over r, and extracting terms proportional to positive powers of r'K This 
reveals potential divergences proportional to r^ and r^, but in both cases the bulk and surface contributions 



*This procedure will not correctly extract divergences of the form log(r/^) . We do not expect to encounter such terms in AL 
theories with z = 2, but in other theories a more careful accounting of the divergences coming from the bulk term is needed. 
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cancel or combine to form a total derivative on the boundary that can be discarded. Thus, the minimal action 
(3.20) is free of r — ^ cxd divergences on-shell. 

Earlier in this section, we mentioned that other sets of surface terms can be used to construct a minimal 
action. This is because the scalar A°'Aa approaches a constant at spatial infinity 

lim A^Aa = -1 . (3.22) 



As a result, the A°'Aa term in (3.4) could be replaced with \/—A°-Aa , as in [24], or some other sufficiently 
well-behaved function. Likewise, the coefficients of the R and F Fab terms could be promoted to functions of 
A°'Aa. Then the surface term in (3.4) would be 



n'^ Jm, 



d^xV^ (^fo{A^A,) + MA" Ac) R + MA'' Ac) F'^'F^b^ , (3.23) 



The calculations in this section place conditions on the on-shell values of the functions fn and their first 
derivatives, and any set of functions that meets those conditions provides a minimal action. Of course, (3.23) is 
by no means the most general set of surface terms that could be used. Terms like R AaAb and F"'cF AaA^, 
where two-derivative functions of the fields are contracted with two-index tensors constructed from Aa, could 
also be included^. Without additional conditions on 61, or some principle that restricts the functional form of 
the action, the surface terms are not uniquely specified. 



3.2 An Extended Action That Allows Independent Variations of a'°' and 



^(o) 



The space of field configurations allowed by a well-defined variational principle must include generic fields that 
differ from solutions of the theory at NNLO. The minimal action derived in the last section enlarges this space, 
to include fields that fall off more slowly than NNLO. Roughly speaking, finiteness of 61 for Hamilton- Jacobi 
variations of the fields implies that 61 vanishes on-shell for field variations that fall off faster than r^, as long 
as 6a = 6(1). The condition on 6a and 6(f) is needed because certain cancellations in the last section relied on 
the variations respecting the constraint on the boundary data. 

It is possible to extend this further, and construct actions such that 61 vanishes on-shell for independent 
variations of the fields that fall off faster than r^ at spatial infinity. Then the space of allowed field configurations 
includes all fields that obey AL boundary conditions, even if they do not admit an asymptotic expansion of 
the form (2.24)-(2.26). This is equivalent to requiring that the action have a finite on-shell response to leading 
order variations of the fields that do not satisfy the constraint on the boundary data. For lack of a better name, 
an action with this property will be referred to as 'extended'. As we will see in the next section, a variational 
principle based on an extended action lets us explore two inequivalent definitions of a boundary stress tensor. 



^Terms that involve four or more derivatives of the fields are not relevant. Compared to the zero- and two-derivative terms, they 
carry additional factors of r^^ which suppress their contributions to the action. 
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An example of an extended action can be constructed by adding another set of surface terms to (3.4) 
2k^ J Mr V 4 '^2 J K'^ Jm, 



^r = WZ2 d''x./^{''R-2A--F^'''F^,-—A^'A^+— jd-'xV-hK (3.24) 



+ \ [ d^X V^ (co + Ci A^Aa + C2R + C3 F'^'^Fab + C4 R^^AaA^ + C5R A^Aa 
K"^ J Mr \ 

+ C6 F'^'Fab A^A, + crF\ F^^ A^Ab 

The on-shell variation of the action takes the same general form (3.5), with the replacement p — t- p + Ap 
and p"" ^ p"" + Aj?*^. The shift in p and p"^ come from varying the new surface terms in (3.24), which changes 
the previous results (3.6) and (3.7) by 

^p^^ = C4 ( ^ /i"^ R^'^ AcAd - R^c A'' A'' - rJ' A^A"" + ^ VcV''{A''A'') + ^ VcV\A''A'') (3.25) 

- i WM^A') - i r" V.VdiA'^A^) 
+ c. ( i k^^RA^A. - RA^A^ - R^^ A.A^ ^ V^W^.) - ^^ V^V.(^^^. 
^c.{U^^F^^F.dA^A.-2F^.F^^A^Ad-F^^F.dA^A^ 
^cr{lk^^F^.F^^A.Ad-F^F^^AdA^-F^F^^A^Ad-F-F^^A.Ad 
Ap" = IciR"'' Ab + 2c5 RA^ + C6 (2 F^'^F^^ A" + 4 F"^ VbiA^A^) - 4 A^A^ VbF^" \ (3.26) 

+ C7I2 F^cF^^ Ab + 2 A^Ab VcF"^ + 2 F"^ Ab V cA" + 2 F"'' A'^ Vc^b 



+ 2 A'^Af, VcF*'^ + 2 F'''^ ^fc Vc^° + 2 F*'^ A'^ Vc^fe ) 



Repeating the calculations of the previous section, with (5a^°' and 5(j)^°^ now treated as independent variations, 
determines the coefficients of the surface terms to be 

'' = ~2l '' = 2i ^^ = "16 "^ = 8 '^=16 "^ = "16 ^^=16-r^- ^^''^^ 

All but one of the coefficients {cj} have been fixed, with cj left undetermined. However, the action turns 
out to be independent of cj when evaluated on-shell. This is because the last two surface terms in (3.24) are 
proportional to each other, up to terms that vanish as r — ^ 00, when evaluated on a solution of the equations 
of motion 

^/^F\ F^"" AaAb l^^o = ^ V^F'^'Fab A' Ac \^^^ . (3.28) 
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This results in a cancellation in the on-shell action when cg is given by (3.27) 



-h [ C6 F-^Fa, A^A, + C7 F\ F^' A^A, 



-h-^F^^F„uA''Ar 



£=0 



16 



(3.29) 



£■=0 



In fact, the contributions that F°'^FabA'^Ac and F"" c F^"^ AaAi, make to Ap'^'' and Ap" satisfy equalities similar 
to (3.28), which leads to on-shell cancellations in 51 that remove terms proportional to cy. Therefore, the last 
surface term in (3.24) will not contribute to any of our calculations, and we are free to set cy = 0. The extended 
action is 



1=^1 d^x^/^ai^R - 2A - If^'^f., - !^^m . ) + 1 



2k^Jm^'^^ ^V^" "^^ 4^ ^^"^ 2'^^^n'K^J9M 



d^xV^K 



(3.30) 



H JdM 



+ / d-'xV-hi - — + — A''Aa-—R+-F'"'Fab + — R'"'AaA, 



16 



-—R A^Aa + — F^'^Fab A" Ac 
16 16 



16 



Like the minimal action (3.20), this action is free of r — )• oo divergences when evaluated on-shell. 



4 Boundary Stress Tensors and Conserved Charges 



Given a variational principle, we can construct a boundary stress tensor and use it to compute the conserved 
charges associated with the asymptotic symmetries of the theory. There are two notions of a boundary stress 
tensor that seem relevant for theories with AL boundary conditions. 

The standard construction due to Brown and York [26] gives the boundary stress tensor as the response of 
the on-shell action to a variation of the boundary metric. This construction can be adapted to non-compact 
spacetimes by working on a compact region Mr: as we did in the previous section, and varying the functional 
/,- with respect to the metric on the boundary M^- 



ab 



5Ir 



-h 6hab 



(4.1) 



£■=0 



Then one can define the charge Q [^] associated with an asymptotic symmetry ^° as the flux of the current Tab C^ 
across a cut C of spatial infinity 



QsYm=Jd^X^cU''Tabe 



(4.2) 



with u" the timelike unit vector normal to C. When the metric is the only field with support at spatial infinity, 
the charges (4.2) are conserved and generate the asymptotic symmetries of the theory. However, the authors of 
[27] showed that this approach must be modified when there are non-vanishing tensor fields at spatial infinity 
in addition to the metric. For a theory with AL boundary conditions, their construction (reviewed below) 
instructs us to replace r°* in (4.2) with an improved boundary stress tensor 



-lab 



t"" + B'^A 



a Ab 



(4.3) 
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where 0° is given by 

1 SI 



-h SAa 



(4.4) 

£=0 



As we will show, (4.3) is the only definition of a boundary stress tensor that is consistent with the minimal 
action defined in section 3.1. Both constructions can be used with the extended action derived in section 3.2, 
and there is an interesting relationship between the charges in that case. 

4.1 The Brown- York Definition Does Not Work For Minimal Actions 

For a variational principle based on the minimal actions defined in 3.1, the leading order variations of hu = —o? 
and A^ = (\> are constrained by (2.22). Since htt and At cannot be varied independently, it is not clear how to 
apply Brown and York's definition of the boundary stress tensor for these actions. 

One might choose to ignore this issue and simply calculate (4.1) for a minimal action, hoping for the best. 
This exercise makes the problem much more concrete. The resulting tensor, when used in (4.2), gives charges 
that are not defined when we take the r — )■ oo limit. To see this in detail, consider the on-shell variation of the 
action 

SI l^^o = \^f^ ^^ ( \ ^"' ^^-^ + ^" ^^- ) • (4-5) 

The coefficients of the field variations are given by 

T°^ = 2{ tt"'' + p"^ ) e'' = TT''+p'' (4.6) 

where p andp"' were defined in (3.6) and (3.7) respectively. Both r** and 0* vanish on-shell (at least to NNLO), 
and the asymptotic expansions for the remaining components are *^ 



r''= (l) rS^ + (-] rd,+ [-) -(i) + ■ ■ • (4-7) 



A' .. fl 



6 



r''=[-) r^^+[^) r^+(-] 4*) + • • • (4-8) 



'^(o) + f^)'^[.)+(i) ^w + --- (4.9) 



The requirement that 81 vanishes on-shell for independent NNLO variations of the fields, as in (3.8), means 
that the leading order terms in all three expansions vanish. The NLO terms in the expansions are relevant 
when we consider the response of the action to the Hamilton- Jacobi variations of the fields (3.14). In general, r^ 
divergences appear in 61 when these terms are non-zero. Since aij can be varied at leading order independently 
of the other fields, finiteness of 51 implies t*2) = 0. But the constraint (2.22), which forces 5a = 5(p, means that 
the NLO terms in r** and 0* do not need to vanish individually. Instead, they satisfy 

l^«(_4e4x)+ 0*^^(2 e2x) = o. (4.10) 



®Note that these expansions are for quantities with raised indices, while most of the other asymptotic expansions throughout the 
paper are for tensors with lower indices. 
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One can check that, generically, tB. and OL. are non-zero. 

The fact that tB. ^ for the minimal action causes problems when we attempt to compute the conserved 
charge associated with the asymptotic symmetry ^°' da = dt using (4.2). With our boundary conditions it is 
natural to take C to be one of the constant t surfaces S^, in which case 

Q[dt]= f d^xV^u'nte ■ (4.11) 

At large r, the factor of ^/a grows like r^, the factor of u* = a^^ falls off as r~^, and the leading behavior of t^ 
is 



Ttt = htthttr'' = l^-j e«>^ r« + . . . . (4.12) 

Combining these factors, we have 

Q[dt]=(^)'' J^d-'xV^e^^Tl^, + ... , (4.13) 

where '. . .' denotes terms of order r^. Thus, attempting to force the Brown- York definition of the boundary 
stress tensor on a minimal action leads to charges that are not defined in the r — )■ oo limit. One might wonder if 
the integral in (4.13) vanishes, leaving only the finite sub-leading part, but this is not the case. For AL solutions 
the integrand is 

e^x r« = e'^ (-^ 7^'°) + {^'x-^^ 'D'x ^^x) • (4.14) 

Even if we restrict our attention to boundary conditions where x is a constant, requiring that the integral of 
7^^°' should vanish is contrary to our original goal of studying asymptotically Lifshitz spacetimes with curved 
spatial sections. 

That the energy has contributions from both the asymptotic part of the metric and the asymptotic part of 
the Proca field has been noted previously [6]. In asymptotically flat or AdS cases, there is enough gauge freedom 
to remove the additional Proca part, whereas for AL spacetimes failure to include the Proca part results in 
a loss of diffeomorphism invariance [6]. In this sense it is not surprising that variations of metric and Proca 
cannot be considered independently. The Br own- York construction is simply the wrong approach for a minimal 
action. 

4.2 The Hollands-Ishibashi-Marolf Boundary Stress Tensor 

It is not surprising that the Brown- York approach fails for minimal actions, since there is a constraint that 
prevents hu from being varied independently of Af. Instead, we must employ the results of Hollands, Ishibashi, 
and Marolf (HIM), who showed in [27] that the construction of conserved charges should be modified when 
other tensor fields besides the metric have support at spatial infinity. 

To define the HIM boundary stress tensor, we first introduce a set of frame fields in the asymptotic region. 
In terms of these fields, the metric and vector are 

Kb = rJAB ^a"^ %^ ^a = ^A ^-a (4-15) 
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for a fixed metric rj^^, and a vector uj^. The HIM boundary stress tensor is then defined as 

1 61 



rpa 



A 



-h Se/ 



(4.16) 

£=0 



Note that there is no factor of '2', since we vary with respect to e^ instead of hab- The on-shell variation of 
the action now takes the form 

51 = [ d^x V^ ( T\ 5e/ + y^ 5ujA ) , (4.17) 

JMr \ J 

with r"4 and V^ related to r"^ and 6" by 

T\ = r-' Vab eh"" + ^'^ o;^ V^ = ^ e/ . (4.18) 

To compute conserved charges, we convert the frame index on T"^ to a spacetime index and then apply the 
usual construction. This gives 

Qhim [^]= J d^x^ u" TaU\ (4.19) 

where T°'^ with two spacetimes indices is 

j,ab ^ j,a^ gb^ ^AB ^ ^ab ^ ffa ^b _ ^4 20) 

This combination of t""^ and 6"" neatly avoids the problems that we encountered trying to force the Brown- York 
construction on a minimal action. With AL boundary conditions, the HIM and Brown- York stress tensors differ 
only in their t-t components 

rptt ^ ^tt j^0t^t ^ ^4 21) 

which is the relevant component when computing the conserved charge associated with the asymptotic symmetry 
generated by dt- The NLO term in the asymptotic expansion of T** vanishes by (4.10), so there is no obstruction 
to defining charges like there was in the Brown- York approach. Expressions for the components of T at NNLO 
are given in appendix D. 

Finiteness of the charges (4.19) follows from (4.10), so it is directly tied to finiteness of 61 for Hamilton- 
Jacobi variations of the fields. In fact, the HIM boundary stress tensor is just the coefficient (up to a constant 
factor) of 5x in the variation of the action. This is perhaps not too surprising, since the kinematics of the 
theory, which in this case includes the constraint (2.22), plays an important role in the analysis of [27]. 

4.3 Brown- York, Revisited 

Although the Brown- York definition of the boundary stress tensor does not work for minimal actions, it can 
be applied to the extended action derived in section 3.2. In that case 61 is finite for independent variations of 
the fields at leading order, which requires r**) = and 0*2) = instead of the condition (4.10). As a result, the 
charges (4.2) are finite as r — )• oo. Furthermore, the fairly restrictive assumptions about the asymptotic form 
of the fields insures that the charges (4.2) are conserved. Expressions for the NNLO terms in the components 
of the HIM and Brown- York stress tensors, for the action (3.30), are given in appendix D. 
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It is interesting to compare the Brown- York and HIM definitions of the energy, since both constructions 
can be appUed to actions hke (3.30). Using (4.20) we have 

Qhim m = Qby m + j d'x^c u'^Aa 9b e' (4.22) 

As before, we wih take C to be one of the spatial surfaces S^, so that n* = a^^. The constraint (2.22) imphes 
u°'Aa — ^ 1 at spatial infinity, so the relation between the two notions of energy becomes 

QuimIC] = Qby[^] + e (4.23) 

with G given by 

9= / d^x^/^et (4.24) 

"'St 

For systems with a thermodynamical interpretation, a better way to think of this is to define a 'chemical 
potential' ^ = —u'^Aa for @ [28]. Then for spacetimes with a horizon at r = r^, regularity of the fields requires 
^('"h) = 0, and the difference in chemical potential between the horizon and spatial infinity is A^ = 0— (— 1) = 1. 
The HIM and Brown- York charges are related by 

Qu,m[C]=Qby[^]+@^^ . (4.25) 

In other words, if we take the Brown- York charge to be the standard internal energy of the system - which is 
a function of - then the HIM charge is the thermodynamic potential that depends instead on the chemical 
potential ^. 



5 The Lifshitz Topological Black Hole 

As an example, we can use the results of the last two sections to compute the action and energy for the 
topological black hole solutions of [7]. Expressed in the coordinate gauge (2.1), the metric and massive vector 
for these solutions takes the form 

g.^dx'^dx'^ = (-Ydr' - l^Xf^irfUirfdt'' + (^Y ^rf fdl^i (5.1) 



r/ \l/ . . . ■ ^£j 

2 



where the functions f±{r) are 



A^dx'^=[^] Uirrdt, (5.2) 



/±(r) = l±^(^)'. (5.3) 



The constant k controls the curvature of the spatial slices, and takes the values ±1 or 0. In those three cases 
the two-dimensional line element dT,^ is given by 

{do'^ + siiPed^'^, k = i 
dd^ + d^d4?, k = Q (5.4) 

^6*2 + sinh^ 6* #2^ k = -I , 
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and the scalar curvature of af-' is 

7^'«) = ^ . (5.5) 

Thus, the /c = solution is just the original Lifshitz spacetime (1.5)-(1.6) written in terms of polar coordinates 
on S(, while the k = —1 solution corresponds to a black hole with a horizon at r^ = ^/v8 and hyperbolic 
spatial slices. The k = \ solution is pathological due to a naked singularity at r = 0, and is therefore only of 
interest as a potential asymptotic completion of some other solution that resolves the singularity. 

In terms of the parameterization of the fields introduced in section 2.2, the topological black hole solutions 
corresponds to AL boundary conditions with x = 0. The NLO and NNLO terms in the asymptotic expansion 
of the fields, which can be taken directly from (5.1) and (5.2), are 

1.2 

(5.6) 
(5.7) 

(5.8) 

These expressions correspond to the results (2.34)-(2.36) obtained at the end of 2.2, with the transverse-traceless 
spatial tensor at NNLO given by Yij = 0. 

First we will consider the minimal action (3.20) for this solution, and calculate the HIM boundary stress 
tensor and conserved charges. The NNLO term in T** is 



-ff 


4 u 


-? = 


64 *J 




a'^' = 


a(^' = 


64 




^" = i 


,/.(^) = 


" 64 ■ 



r« = ^^ in^^^f = ^ , (5.9) 

(4) 128 ^ '^ 32^ ' ^ ^ 
so the energy obtained from the HIM definition of the conserved charges is 

QmM[dt] = f d^x^/^T'i) (5.10) 



evoUJ:,) ^^ 



2,21k? ' 

where Volfc(Si) is the (dimensionless) volume of the surface Sf with metric (5.4). The energy and on-shell 
action are both zero for the A; = solution, as expected, so let us focus on the case k = —1. This solution has 
a horizon at Th = ^/\/8, and regularity of the metric implies a periodicity /3 = T"^ = Ani in Euclidean time. If 
we take the entropy S to be one-quarter of the horizon area in Planck units, then we find 

rs = £2voU(Si)^. (5.11) 

Evaluating the action (3.20) on the Euclidean section of the solution gives 

I, = -£2voU(S0^^, (5.12) 

which satisfies the expected relation between the action, energy, temperature, and entropy 

h = /3{QniM[dt]-TS) . (5.13) 

23 



It is important to point out that both the energy (5.10) and the on-sheh action (5.12) are sensitive to the choice 
of surface terms in the minimal action, but the result (5.13) is not. 

Now consider the extended action (3.30), which contains additional surface terms compared to the minimal 
action. These terms shift the value of the energy obtained via the HIM construction, which is now 

QniMidt] = fVohi^t) ^^ . (5.14) 

However, the action evaluated on the Euchdean section is shifted by the same amount relative to (5.12), so the 
relation (5.13) is still satisfied. For the k = solution the action vanishes, as before, while for the k = —1 it is 



7_ 
32Tk2 



-fYol^i{j:t)i;:TTZ2 ■ (5.15) 



Of course, the action (3.30) also lets us calculate the Brown- York stress tensor, which is given by 

rft = -1^ (7^(«))2 + -^ Di!D*7^(°' - -^ </.'*' . (5.16) 

<*' 128^2 ^ ^ 16^2 * £«;2 ^ ^ ^ 

It is interesting that, unlike the HIM boundary stress tensor, this depends on the part of the solution at NNLO 
that is not fixed by the kinematics of the theory. The charge associated with dt obtained from this stress tensor 
is 

1.2 

QBY[5i]=^2Volfc[Si]— . (5.17) 

To relate this to the HIM energy, we must also evaluate the quantity G in (4.24). The NNLO term in the 
asymptotic expansion of 0* is 

^(4) = ^ (^'°f + ^ ^«2)^7^'°) - -^ r^ , (5.18) 

which also depends on the dynamical part of the solution. Then is 

G = -|^ d'xV^el, = -fVoU^t) ^^ ■ (5.19) 

The chemical potential for G is ^ = —uP'Aa, which takes the form 

* = --<^ = -TTT (5.20) 

a f-{r) 

for the solution (5.1)-(5.2). When k = —1 this is regular on the Euclidean section, vanishing at r^ and 
approaching — 1 as r — >■ oo. Thus, the Brown- York and HIM charges satisfy 

Qu,M[dt]=QBY[dt] + @A^ , (5.21) 

where A^ = ^{r^) — \I'(oo) is the difference in chemical potential between the horizon and spatial infinity. 



6 Discussion 

There are several ways in which our analysis can be improved. The most obvious extension is a treatment 
of asymptotically Lifshitz spacetimes that does not rely on the assumptions (2.12)-(2.13). Relaxing these 
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conditions is challenging, since the 3—^2 + 1 split becomes much more complicated, but this is just a technical 
difficulty and we do not anticipate any real difficulties. A broader treatment of AL spacetimes would have 
a significant impact on many of our results. First, there will be a larger class of field variations consistent 
with more general boundary conditions. Requiring 61 to have the appropriate properties for the full class of 
field variations will give additional conditions on the surface terms, and this may partially or entirely resolve 
the ambiguities encountered in section 3. A second consequence of more general boundary conditions is that 
the Brown- York charges, if they can be defined, will no longer be conserved. The Hollands-Ishibashi-Marolf 
construction [27] should work perfectly well in that case, but the charges will not admit an interpretation along 
the lines of (5.21). 

Another relevant question is how our results extend to other values of the critical exponent z. The main 
difficulty here is that the character of the asymptotic expansions depends on z. For example, NNLO terms and 
NLO-squared terms enter the asymptotic expansions with different powers of r when 1 < z < 2, but with the 
same power of r when z = 2. This issue affects every aspect of the analysis in sections 2 and 3, making the 
prospect of a general result for arbitrary values of z seem unlikely. Instead, we expect to find actions that apply 
for distinct ranges of the critical exponent where functions of the fields have qualitatively similar asymptotic 
expansions. It would be interesting to see if there is a suitable action for some range of z that includes z = 1, 
since this corresponds to asymptotically AdS4 spacetimes where the massive vector does not have support at 
spatial infinity. The surface terms in the action are unique and well-understood in that case [18, 19]. 

A further avenue of study would be to understand the relationship between our results and the stability of 
AL spacetimes. An initial value analysis recently showed that a generic normalizable state in an AL spacetime 
will evolve in such a way to violate Lifshitz asymptotics in finite time [6]. Whether or not the appropriate 
counterterm action can shed further light on this subject remains to be seen. 

Finally, we have focused on the action as it relates to the gravitational theory, emphasizing the criteria for a 
well-defined variational problem and the definition and application of various boundary stress tensors. Despite 
the fact that we were originally motivated by a duality between gravitational and condensed matter theories, 
we have not discussed our results as they might apply to the CM dual. Once we have obtained a sufficiently 
general definition of AL boundary conditions, it will be interesting to see how the ambiguities in the action - 
if they still remain - relate to the structure of 1-, 2-, and higher n-point functions of operators in the dual CM 
theory. 

Acknowledgments 

This work was supported in part by the Natural Sciences and Engineering Research Council of Canada, and by 
Loyola University Chicago's Summer Research Stipend Program. Many of the calculations in this paper were 
checked using xAct [29], a tensor computer algebra package for Mathematica. 

Note: As this paper was being prepared for submission, the papers [30, 31] appeared on the arXiv. They take 
a systematic approach to holographic renormalization, along the lines of [32, 22], for general values of z and 
with a different definition of asymptotically Lifshitz boundary conditions than the ones we present in section 2. 
Our results seem to agree with theirs where they overlap. 

25 



A The 4—^3 + 1 Decomposition 

In section 2.1 we carry out a 3 + 1 split of the fields and equations of motions, working in the coordinate 
gauge (2.1). Projections orthogonal to the surface Mj. are obtained by contracting with the unit normal vector 
n^ = {£/r)5'^r, and projections tangent to Mr with P^^ = dx^/dx"". The metric (first fundamental form) on 
Mr is given by the projection of the spacetime metric 

hab = Pa^Pb^g^.. , (A.l) 

and the extrinsic curvature (second fundamental form) is proportional to the Lie derivative of hab along the 
normal 

Kab = Pa^Pb'' 'V^n, = ^ £nhab ■ (A.2) 

The projection Pa^ naturally defines a covariant derivative Va that acts on tensors tangent to Mr and is 
compatible with the metric hab- It is obtained from the complete projection (tangent to Mr) of '^V^ acting on 
the tensor. For instance, given a vector U^ such that U^n>^ = 0, the derivative is 

VaUb:=Pa^Pb''^V^U,. (A.3) 

For an arbitrary 4-vector V^ with components Vn = n^V^ and Va = Pa^ V^ we have 

n^n'' V^K = -fnK (A.4) 

PiTn^ V^K = £nVb - Va K\ (A.5) 

Pa^Pb" V^K = ^aVb + Vn Kab (A.6) 

^V^y^ = i:„T4 + VnK + VaV . (A.7) 

Likewise, for the anti-symmetric tensor F^i, with components given by (2.2)-(2.3), the relevant projections are 

n" V^F^, = - VaB'' (A.8) 

n" 'V^F^, = £r,Bb + BbK-2BaK\ + V'Fab . (A.9) 

These results are specific to the coordinate gauge (2.1). In a generic coordinate system they would acquire 
additional terms involving the vector Pa^ n'^ ^^un^ and its derivatives. 

To derive the 3 + 1 equations of motion we also need the projections of the four-dimensional Ricci tensor 
normal and tangent to Mr- They can be written in terms of the extrinsic curvature Kab and the intrinsic Ricci 
tensor Rat as 

n^n^X-= -£nK-K''^Kab (A.IO) 

Pa" n'^ ^R^.u = y^Kab - VaK (A.ll) 

Pa^Pb'' ''Rf.. = Rab - £nKab -KKab + 2 Ka" K^c ■ (A. 12) 

Using these projections to rewrite the trace of the Ricci tensor gives the four-dimensional Ricci scalar in terms 
of the intrinsic and extrinsic curvature 

^R = R-K^ -K^^Kab-2£nK . (A.13) 
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B The 3—^2 + 1 Decomposition 

The 2 + 1 split in section 2.1 involves projections orthogonal to S^ with the (forward-pointing) timelike unit 
vector n" = a~^ d'^t, and tangent to S^ with Pj" = dx'^/dx^. The results are greatly simplified by the restrictions 
placed on the field in (2.12)-(2.13). In particular, the extrinsic curvature of Sj C Mr vanishes 

% = -Pi^'Pj^ VaUb = -^ (dtcTij - D,pj - Djl3^) = , (B.l) 

where Di is the two-dimensional covariant derivative that acts on tensors tangent to Sj. On the other hand, 
the acceleration vector u^ ^bUa is non-zero if the lapse varies over S^ 

Pi''(u^VbUa) = ^Dia. (B.2) 

The projections of the three-dimensional Ricci tensor are 

u'^v!' Rab= -DiD'a (B.3) 

a 

P^u'Rab = (B.4) 



1 

a 



P^ P/ Rab = Tlij - - DiDja , (B.5) 



where TZij is the two-dimensional Ricci tensor for the metric aij. The three-dimensional Ricci scalar is 

R = n-- DiD'a . (B.6) 

a 

The remaining projections of fields and covariant derivatives are straight-forward. 

C Transverse-Traceless Tensors Constructed from the Boundary Data 

The integral of the Ricci scalar is a topological invariant in two dimensions. As a result, the Einstein tensor 
vanishes identically in two dimensions 



nij--aijn = o (c.i) 



for any metric dij. Now consider a metric a^j + ejij, where ^ij is some well-behaved but otherwise arbitrary 
tensor, and e is a small parameter. Expanding the Einstein tensor for this metric in powers of e gives 

G,, [a + ej]= Gij [^] + | ( D'^Dakj + D^D.j^,k ~ D^D^jij - D,Dj^\ - IZ 7,,- (C.2) 

where D^ is the covariant derivative compatible with (Tjj, and indices are lowered and raised using aij and its 
inverse. The Ricci tensor vanishes identically in two-dimensions, so the left-hand side of (C.2) and the first 
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term on the right-hand side are both zero. And since e is a continuous parameter, the remaining terms in the 
Taylor expansion must vanish order- by-order. For the 0{e) term, this means 



= D'^Dakj + D'^Dj^ik - D^DkJ^J - DiDjj\ - n ^ij (C.3) 

In other words, the hnearized Einstein tensor vanishes for a smah perturbation of the metric. 

The result (C.3) is useful if we regard it as an identity that holds in two dimensions for a metric aij and 
an arbitrary symmetric tensor 7jj . If we consider different tensors 'jij constructed from derivatives of a scalar 
function ip, we obtain a set of identities that are useful when solving the equations of motion at NNLO in 
section 2.2. First, let 7jj be 

7i,- = D^^pDjiP . (C.4) 

Substituting this in (C.3), commuting covariant derivatives, and applying (C.l), we obtain 

= 2 D^Djip D'^Dkip - 2 DiD'^Tp DjU^Tp + aij D'^D^ip DkDiip - aij D^Dkip D^Diip . (C.5) 

At first glance this combination of terms is not obviously zero, but the result can be confirmed by direct 
calculation. Another useful identity comes from setting 

jij = e""^ DiijDjij , (C.6) 

with c a constant. Then (C.3), supplemented with the identity (C.5), gives 

= DiTp Dj^j D^Dk^ + DiDji) D'^Tp Duip - DiDkip Djip D^i) - DjDui) Diip D'^^j (C.7) 

- aij D^Dkip D^^ DiiP + aij D^D^^ Dk^ DuP . 

It is easy to check that the expressions on the right-hand sides of (C.5) and (C.7) are traceless. One can also 
check that they are transverse, which requires commuting covariant derivatives and multiple applications of 

(C.l). 

There are two useful applications of the identities (C.5) and (C.7) in section 2.2. First, solving the equations 
of motion for the NNLO term in the spatial metric is complicated, and applying identities like (C.l) at different 
stages in the calculation leads to expressions for aPPj which differ by linear combinations (C.5) and (C.7). 
Second, these identities insure that the transverse-traceless tensor Yij that appears at NNLO in the spatial 
metric does not depend on the boundary data. Consider a generic two-index symmetric tensor constructed 
from four-derivative combinations of aij and tp, with at least one derivative acting on each factor of ip. Thus, 
terms like DiipDjipD^Di^tp or aijTZD^Dj^ip may appear, but not ipDiDjipD^Dj^ip or aijip D^Dj^TZ. The only 
transverse and traceless tensors of this sort are given by the terms on the right-hand sides of (C.5) and (C.7), 
which vanish. Therefore, there are no transverse-traceless four-derivative functions of the boundary data that 
can appear in af^ . 
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D NNLO Terms in the Boundary Stress Tensors 

We present here expressions for the components of the boundary stress-tensors for the minimal action (3.20) 
and the extended action (3.30). 

HIM boundary stress tensor 

The NNLO terms in the components of the HIM stress tensor for the minimal action (3.20) are: 

- -^ V'xI^kX ^'x^ix - Y ^"^'^ '^kX 'Dix - Y ^'^'^ ^'^^'^ 



£3 .. ^3 ^3 n3 

16 



4 7^(°) !D\ 'Dkx - -^ n'"' '^\ T^'x - -^ m^V'Vkx + -^ 4 '^'"Dki)''^ix 



+ -^ 4 D'^D'x 2)fcx 2)/x - ^ 4 ^'2?fcx i)'x i)a + Y ^'^ ^'^ ^'^ ^'^■^ 



'(0) 

For the extended action (3.30), the additional surface terms shift the NNLO terms in T . In that case they 
are given by: 

Tl, = e-^^ ( ^ (7^(°))2 + i! !D'=2)fc7^(°' + ^ 7^(°) D'^ID.x + ^ I)fc7^(°) D\ - -^ 7^(°) D'^xI^fcX (D.3) 
+ 1^- !D^D,x I^'Da + ^ I^'2);iX I^'x 2^a + -^ ^"^ 2)'2),D,x + ^ O^'^IJfcD'lJ.x 
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'^ = iS '^ ^^'"'^^ " ^ ^^^'^'°' + ^ ^(0) 25'=2)fc7^(°' - ^ 7^(°) O^^O^^x (D.4) 



+ ^ 4 n^'' 2)\ 25fcX - ^ T^'"' 2)^x 25^X - ^ 2)*a?-'D'=DfcX + ^ 4 a^^'a^^D^D^x 



- ^ (d^x O^^'a^'^D.x + "^'X 'D"D'"Dkx) + ^ 4 '^kX 'D^'D'x 

- -^ !D*D^x ^'■DkX + ^ 4 2)'2),x 'D^'DiX 

+ £3 D^Di^ dA^;^ o^.x - -^ (d'x ^^^'DfcX + 'D^x '^''Dkx) T^kX 

+ -^ 4 I^'^D'x 25fcX I)a - ^ 4 ^'2),.x D'x 1>iX + Y ^'X 2)^'X O^'^x l^^A^X 

-^4D^x2)fcx2)'x2)a. 

Brown- York boundary stress tensor 

The Brown- York stress tensor can also be defined for the extended action (3.30). In that case the component 
r*4) is the same as (D.4), and the t-t component is given by 



+ ^ 7^(°) D'^xI^fcX + ^ -D'-^kX 25'D/X + ^ i^'^DfcX 'D'X '^iX + ^ i^'^x I^'l^z^^fcX 
+ 4r D'DfcD'Da - 4r 2)'xI)fcX I^'x^^a - ^ I^'^^'x I^fcX I^a - ^ I^'^^'x DfcO^/x ) ■ 



We also require the NNLO term in the expansion of O"^. With AL boundary conditions the component 0* 
vanishes, and 6\^^^ is 

^''^ = "'''' ( ^ ^^''^' " I '^''' + ^ a3'^Dfc7^(°) + -^ Dfcx 2)'7^(°) + -^ 7^(«) v^ '^kx (d.6) 

--^D^D'xI^fcl^a') ■ 
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